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1. Introduction 

The discovery of cosmic acceleratiorM^ has stimulated an intensive study of mod- 
els of dark energy^ responsible for the origin of this phenomenon. Dark energy is 
characterized by a negative pressure whose relation to the energy density w = p/e 
C^ ' is less than —1/3. Moreover, if this relation happens to be less than —1 such kind 

of dark energy is called "phantom" dark energy^. The breakdown of the condition 
w > —1 implies the so called super- acceleration of cosmological evolution, which in 
some models culminates approaching a new type of cosmological singularity called 
Big Rip^. While the cosmological constant {w — —1) is still a possible candidate 
for the role of dark energy, there are observations giving indications in favor of 
the models where the equation of state parameter w not only changes with time, 
but could become less than — 1 no^. There are also observations '^ indicating the 
phenomenon of the phantom divide line crossing '^, when the parameter w changes 
between the phantom and non-phantom phases. A standard way to introduce the 
phantom energy is to consider a scalar field with the negative sign of the kinetic 
energy term, usually called "phantom" scalar field. The simplest way of describing 
the phenomenon of the phantom divide line crossing is to introduce two scalar fields 

Q 

- phantom and non-phantom as in the so called qui ntom mo dels . Some other mod- 
els are based on the use of the PT symmetry ideas l ^^ i Note that the phantom 
divide line crossing can occur also in models based on one scalar field if some special 
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potential are chosen Lt^ or if the non-minimal coupling between a scalar field and 
gravity is included L!^. 

Naturally, consideration of two scalar fields open wider opportunities for de- 
scription of any prescribed cosmological evolution, than working with a single field 

n 

Recently, another class of cosmological dark energy models was put forward 
- the so called chameleon cosmology models ^^. In this models the scalar field 
interacts with some kind of matter, behaving as perfect fluid. The Lagrangian of the 
model contains a term, where the matter interacts with the effective metric (physical 
metric, multiplied by a conformal factor depending on the scalar field). The energy 
density of the perfect fiuid is proportional to the trace of its energy-momentum 
tensor, multiplied by this conformal factor. On the other hand the effective mass of 
the scalar field becomes variable, depending on the energy density of matter, which 
changes in the process of the cosmological evolution. 

Lately, an alternative version of the Chameleon cosmology was proposed in paper 
Li^. Here the term, describing the interaction between the chameleon scalar field 
and matter has a very general multiplicative form and is equal to some arbitrary 
function of the scalar field multiplied by the Lagrangian of matter, which on the 
Friedmann- Robertson- Walker background coincides with the pressure of the matter. 
Thus, one has two arbitrary functions of the scalar field to play with - the standard 
potential and the potential-like function describing the interaction of matter with 
the chameleon field. Choosing properly these two functions of the scalar field one 
can describe different kinds of cosmological evolution, including those where the 
universe undergoes the phantom divide line crossing. The authors of paper ^^ use 
the numerical integration of the corresponding equations of motion and compare 
the results with the data coming from time evolution of the cosmological redshift 
of distant sources. 

We remark that having two functions of a scalar field one gets a sufficient freedom 
to construct exactly solvable solutions of the cosmological equations, giving some 
prescribed cosmological evolution. Just like choosing exponential potentials one can 
find exact solutions desribing power-law expansion of the universe (see e.g. ^^^, 
here one can try to construct two potential-like functions providing the existence of 
exact cosmological solutions, undergoing the phantom divide line crossing. In fact, 
this is the task of the present paper. 

We consider the flat Friedmann model with the metric 

ds^ ^dt^ -a'^{t)df. (1) 

The Hubble parameter behaves as 

^ ' a{t) t{tR -ty ^ ^ 

where a is some positive parameter of the order 1 and "dot" means as usual the 
derivative with respect to the cosmic time t. Such a cosmological evolution was 



May 12, 2010 0:14 WSPC/INSTRUCTION FILE chameleon 



Chameleon cosmology model describing the phantom, divide line crossing 3 

considered in papers 1^ 1 I in the the context of two-field models and was called 
"Bang-to-Rip" evolution. Here the universe begins its evolution from the Big Bang 
singularity when the Hubble parameter behaves as H{t) ^ j and ends in the Big 
Rip singularity at the moment t = tn. At the moment t — ^ the phantom divide 
line is crossed by the universe. Indeed, the phantom divide line crossing means 
that the relation between the pressure and the energy density becomes equal to —1. 
That means that in the energy conservation law i ~ —3H{e+p) the right-hand side 
vanishes, and so does the time derivative of the total energy density in its left-hand 
side. The energy density is proportinal in the flat Friedmann model to the Hubble 
parameter squared. Hence, the time derivative of the Hubble parameter H should 
vanish too. The time derivative of the Hubble parameter from Eq. ([2|) is 

. 2atR{2t-tR) 

""= tHtn-ty (^) 

and vanishes at t — tfl/2. 

We shall consider the minimally coupled to gravity scalar field 4> with the po- 
tential V{(j)) whose interaction with the perfect fluid is descr ibe d by the term in the 
Lagrangian, which on the Friedmann background looks like^!^ 

Lscalar+matter = lPf{(l>): (4) 

where the coefficient 7 relates the energy density p and the pressure p of matter: 

P ^ IP, (5) 

and f{(j)) is some function of the scalar field cj). In the next section we shall find the 
form of the functions V{(f)) and /(</>), and the exact solution for a{t) and 0(i) match- 
ing the evolution ([2]). In the third section we discuss the bounds on the parameters 
of the model. The last section contains some concluding remarks. 

2. An exactly solvable Chameleon field cosmological model 

We shall fix the fundamental constants in such a way to give to the Friedmann 
equation a particularly simple form: 

H^ = e, (6) 

where e is the total energy density of the scalar field and matter. On the fiat 
Friedmann background this total energy density is 

e^^ + V{^) + pf{^). (7) 

The Klein-Gordon equation for the scalar field is 

^ + 3i/0+l^'(0)+7P/'(0)=O, (8) 

where "prime" stays for the derivative with respect to <j). The total energy density 
e satisfies the energy conservation law 

e + 3iJ(e + P) = 0, (9) 
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where P is the total pressure of the matter and of the scalar field which is equal to 

p = 7P+Y-n0)- (10) 

Substituting Eqs. ^ and ^ into Eq. (O we obtain 

(A(^ + 3i/0 +V' + pf) = -3(1 + ^)Hpf - pf. (11) 

Using Eq. ([8]) we can simplify the left-hand side of Eq. (|lip coming to 

(l-7)p/ + 3(l + 7)i^P/ + P/ = 0. (12) 

Integrating Eq. (|12p one obtains for the energy density of the matter 

^^*^ " /i-^(</)(i))a3(i+^)(i)' ^^^^ 

where po is an integration constant. 

Now the Friedmann and Klein- Gordon equations can be rewritten as 

4, + 3H4, + V' + j^^e^ = 0. (15) 

We have chosen the cosmological evolution described by the Hubble variable ([2]) . 
Correspondingly the cosmological radius a{t) is 

where oq is a positive constant. Now we can rewrite the matter energy density p(t) 
as 

Af(t^-t)3"(i+7) 
/'W ^ ;i-7t3a(i+7) ' (1^) 

where we have introduced a new constant 

M=— ^. (18) 

3(1+7) ^ ' 

"0 

The term pf in the right-hand side of the Friedmann equation (J14p now looks as 

MpjtR - t)3"(i+^) 
^/ = -^^T+T) • (19) 

Now, the left-hand side of Eq. ([T4| behaves as ~ l/t^{tR — t)"^. If we would 
like the contribution of the kinetic term of the scalar field 0^/2 to have the same 
structure we should require that the time dependence of the scalar field is 

0(t) = <^oln-^, (20) 

IR -t 
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where 0o is a constant to be determined. Then, we shall look for the time dependence 
of the potential V in the form 

We shall try to find time dependence of the function / from the relation 

r = (/o + /it)t'"('+^)-'(ii? - i)-3"(i+^)-2. (22) 

Substituting Eqs. ^, ^, ([21]) and ([22]) into Eq. ^ we obtain the following 
constraints on the constants : 

ahl^Mfo + ^ + Vo, (23) 

Vi + A//i = 0. (24) 

Now, let us consider the Klein-Gordon equation ([T5|) . The second time derivative 
of the scalar field is 

7 __ _(h 4>o _ 4>otR{2t - in) 

^~ t^ {tR-tr ^ t^{tR-tY ■ 
The "friction" term is 

3g'^- ,2(,^_,)2 - (26) 

The derivative of the potential is 

and the term, including the interaction with matter is 

1 

(j^otRt^JtR - tf 

+ (4/o + /iifl,(3a(l + 7) - 1) + S/it^). (28) 

Now, substituting the expressions (IT5|) - (P5)) into Eq. (IT5|) we obtain the following 
constraints: 

-4>ltl + ?,4>ltla - 2Vb + M/o(3a(l + 7) - 2) = 0, (29) 



(25) 



^f'P = ^..2U „,^2 X (/oifl(3a(l + 7) - 2) 



20^<^ + 4Vo - YxtR + 4Af /o + M/it;^(3a(l + 7) - 1) - 0, (30) 

Wx + 3M/i = 0. (31) 

The last condition ([5T|) coincides with ([Ml) . 

Now we can solve the constraints and find the relations between the parameters. 
We shall choose as free parameters the amplitude of the scalar field (/)o and the 
"quantity of matter " M > 0. From Eqs. ^ and ^ we find 
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Then from Eq. dM]) fonows 



^' ' 'imr^y ^^^^ 



From Eqs. ^ and ^ we find 



^° 3M(l + 7) ^^^^ 



and 



,r _ t|(3a(3a(l+7)-2) + 3</>g(l-7)) .„. 

We should also now express the time parameter t in terms of the scalar field 0. 
It follows immediately from Eq. ([20]) that 

t= —, r- (36) 



1 + -P(-Il) 



Substituting the parameters Vo,Vi,/o and /i from Eqs. (|32l) -(V0) and the ex- 
pression for the time parameter t from Eq. (|36|) into Eqs. (f2T|) and ([22]) we find the 
explicit expressions for the potential V{4i) and the function f{4i): 

8cosh*^ / ^ ^ oi \ 

n'^)- 3(1 + ^^° (^6a^(l + 7) + 3</>^(l-7)+4«tanh^j, (37) 

1 

16cosh''^expf3a(l + 7)^') / j, 



3. Parameters of the model 

The parameters a and tn. are given by the cosmological dynamics of the model 
(HI). The parameter 4>q is free. Without loosing generality we can choose it positive. 
Correspondingly the value of the scalar field </> (see Eq. ([20l) ) runs from —infty 
at the Big Bang beginning of the cosmological evolution to +c)o at its end at the 
moment t = tn when it encounters the Big Rip singularity. It crosses the zero value 
at i = tfi/2, i.e. at the moment of the phantom divide line crossing. Now, look- 
ing at the expression for the function /((/>) describing the interaction between the 
chameleon scalar field and the matter, we see that it is well defined only if the 
expression, which should be elevated into the power I/7 is always positive, or if the 
parameter 7 satisfies some additional restrictions. Thus, we can consider two cases. 



The case I. The expression in the bracket in Eq. (|38p is always positive. It is 
possible if and only if 

(39) 
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7 < -1. (40) 



Then if the parameter a > :t and 



3(1-7) 
the potential V{(l>) is always negative. If a > ^ and 



^ Ua-ea^ilpl (41) 




2« ^^ ^ /4a-6a2(l + 7) ^^^^ 



3(1-7) 
the potential V{<j>) changes sign at 

= 200 arctanh5^^!ii±^i±M(l^. (43) 

4a 

If a < i the potential is always negative. 

The case II. If at least one of two inequalities ([39l) , pO]) is broken the expres- 
sion for /''' in Eq. (|38p cannot be always nonnegative. Hence, we should impose a 
following condition on the factor 7: 

7= , (44) 

n 

where m and n are integers. In this case the expression for / is well defined. The 
sign of the potential depends on the interplay of three parameters (pQ, 7 and a. 

4. Conclusion 

We have constructed an exact solution for a particular Chameleon cosmological 
model of the type considered in paper ^^. The universe in this solution begins its 
evolution from Big Bang singularity, undergoes a phantom divide line crossing and 
ends in the Big Rip singularity. The two potential-like functions of the chameleon 
scalar field have a rather simple analytic form. Note that this form is simpler than 
the pote ntial functions in two-scalar model, providing the same cosmological evolu- 
tion LSJJI. Thus, we hope that our version of the Chameleon model can be useful for 
the analysis of the phantom divide line crossing phenomenon. 
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